Set up environment.

Fraxinus americana

speciesAb_ <- 'Frax_Amer'
load(paste0('./Models/model_outputs/', speciesAb_, '_GCM', gcm, 
                            '_PC', pc, '.rData'))
plot(rangeMap, border = 'blue', main = substitute(paste('Maxent output, ', 
                                                        italic('Fraxinus americana '),
                                                        'occurrences')), 
     sub = paste0('GCM = ', gcm))
plot(envMap, add = TRUE)
plot(rangeMap, border = 'blue', add = TRUE)
map("state", add = TRUE)
map("world", add = TRUE)
points(records$longitude, records$latitude, pch = 16, cex = 0.6, col = 'red')

Fraxinus anomala

speciesAb_ <- 'Frax_Anom'
load(paste0('./Models/model_outputs/', speciesAb_, '_GCM', gcm, 
                            '_PC', pc, '.rData'))
plot(rangeMap, border = 'blue', main = substitute(paste('Maxent output, ', 
                                                        italic('Fraxinus anomala '),
                                                        'occurrences')),
     sub = paste0('GCM = ', gcm))
plot(envMap, add = TRUE)
plot(rangeMap, border = 'blue', add = TRUE)
map("state", add = TRUE)
map("world", add = TRUE)
points(records$longitude, records$latitude, pch = 16, cex = 0.6, col = 'red')

Fraxinus caroliniana

speciesAb_ <- 'Frax_Caro'
load(paste0('./Models/model_outputs/', speciesAb_, '_GCM', gcm, 
                            '_PC', pc, '.rData'))
plot(rangeMap, border = 'blue', main = substitute(paste('Maxent output, ', 
                                                        italic('Fraxinus caroliniana '),
                                                        'occurrences')),
     sub = paste0('GCM = ', gcm))
plot(envMap, add = TRUE)
plot(rangeMap, border = 'blue', add = TRUE)
map("state", add = TRUE)
map("world", add = TRUE)
points(records$longitude, records$latitude, pch = 16, cex = 0.6, col = 'red')

Fraxinus cuspidata

** Check values in E Texas. plot(raster) extract(raster, click(n=1))

speciesAb_ <- 'Frax_Cusp'
load(paste0('./Models/model_outputs/', speciesAb_, '_GCM', gcm, 
                            '_PC', pc, '.rData'))
plot(rangeMap, border = 'blue', main = substitute(paste('Maxent output, ', 
                                                        italic('Fraxinus cuspidata '),
                                                        'occurrences')),
     sub = paste0('GCM = ', gcm))
plot(envMap, add = TRUE)
plot(rangeMap, border = 'blue', add = TRUE)
map("state", add = TRUE)
map("world", add = TRUE)
points(records$longitude, records$latitude, pch = 16, cex = 0.6, col = 'red')

Fraxinus greggii

speciesAb_ <- 'Frax_Greg'
load(paste0('./Models/model_outputs/', speciesAb_, '_GCM', gcm, 
                            '_PC', pc, '.rData'))
plot(rangeMap, border = 'blue', main = substitute(paste('Maxent output, ', 
                                                        italic('Fraxinus greggii '),
                                                        'occurrences')),
     sub = paste0('GCM = ', gcm))
plot(envMap, add = TRUE)
plot(rangeMap, border = 'blue', add = TRUE)
map("state", add = TRUE)
map("world", add = TRUE)
points(records$longitude, records$latitude, pch = 16, cex = 0.6, col = 'red')

Fraxinus latifolia

speciesAb_ <- 'Frax_Lati'
load(paste0('./Models/model_outputs/', speciesAb_, '_GCM', gcm, 
                            '_PC', pc, '.rData'))
plot(rangeMap, border = 'blue', main = substitute(paste('Maxent output, ', 
                                                        italic('Fraxinus latifolia '),
                                                        'occurrences')),
     sub = paste0('GCM = ', gcm))
plot(envMap, add = TRUE)
plot(rangeMap, border = 'blue', add = TRUE)
map("state", add = TRUE)
map("world", add = TRUE)
points(records$longitude, records$latitude, pch = 16, cex = 0.6, col = 'red')

Fraxinus nigra

speciesAb_ <- 'Frax_Nigr'
load(paste0('./Models/model_outputs/', speciesAb_, '_GCM', gcm, 
                            '_PC', pc, '.rData'))
plot(rangeMap, border = 'blue', main = substitute(paste('Maxent output, ', 
                                                        italic('Fraxinus nigra '),
                                                        'occurrences')),
     sub = paste0('GCM = ', gcm))
plot(envMap, add = TRUE)
plot(rangeMap, border = 'blue', add = TRUE)
map("state", add = TRUE)
map("world", add = TRUE)
points(records$longitude, records$latitude, pch = 16, cex = 0.6, col = 'red')

Fraxinus pennsylvanica

speciesAb_ <- 'Frax_Penn'
load(paste0('./Models/model_outputs/', speciesAb_, '_GCM', gcm, 
                            '_PC', pc, '.rData'))
plot(rangeMap, border = 'blue', main = substitute(paste('Maxent output, ', 
                                                        italic('Fraxinus pennsylvanica '),
                                                        'occurrences')),
     sub = paste0('GCM = ', gcm))
plot(envMap, add = TRUE)
plot(rangeMap, border = 'blue', add = TRUE)
map("state", add = TRUE)
map("world", add = TRUE)
points(records$longitude, records$latitude, pch = 16, cex = 0.6, col = 'red')

Fraxinus profunda

speciesAb_ <- 'Frax_Prof'
load(paste0('./Models/model_outputs/', speciesAb_, '_GCM', gcm, 
                            '_PC', pc, '.rData'))
plot(rangeMap, border = 'blue', main = substitute(paste('Maxent output, ', 
                                                        italic('Fraxinus profunda '),
                                                        'occurrences')),
     sub = paste0('GCM = ', gcm))
plot(envMap, add = TRUE)
plot(rangeMap, border = 'blue', add = TRUE)
map("state", add = TRUE)
map("world", add = TRUE)
points(records$longitude, records$latitude, pch = 16, cex = 0.6, col = 'red')

Fraxinus quadrangulata

speciesAb_ <- 'Frax_Quad'
load(paste0('./Models/model_outputs/', speciesAb_, '_GCM', gcm, 
                            '_PC', pc, '.rData'))
plot(rangeMap, border = 'blue', main = substitute(paste('Maxent output, ', 
                                                        italic('Fraxinus quadrangulata '),
                                                        'occurrences')),
     sub = paste0('GCM = ', gcm))
plot(envMap, add = TRUE)
plot(rangeMap, border = 'blue', add = TRUE)
map("state", add = TRUE)
map("world", add = TRUE)
points(records$longitude, records$latitude, pch = 16, cex = 0.6, col = 'red')

Now, for Lorenz ecbilt gcm.

Fraxinus americana

speciesAb_ <- 'Frax_Amer'
load(paste0('./Models/model_outputs/', speciesAb_, '_GCM', gcm, 
                            '_PC', pc, '.rData'))
plot(rangeMap, border = 'blue', main = substitute(paste('Maxent output, ', 
                                                        italic('Fraxinus americana '),
                                                        'occurrences')),
     sub = paste0('GCM = ', gcm))
plot(envMap, add = TRUE)
plot(rangeMap, border = 'blue', add = TRUE)
map("state", add = TRUE)
map("world", add = TRUE)
points(records$longitude, records$latitude, pch = 16, cex = 0.6, col = 'red')

Fraxinus anomala

speciesAb_ <- 'Frax_Anom'
load(paste0('./Models/model_outputs/', speciesAb_, '_GCM', gcm, 
                            '_PC', pc, '.rData'))
plot(rangeMap, border = 'blue', main = substitute(paste('Maxent output, ', 
                                                        italic('Fraxinus anomala '),
                                                        'occurrences')),
     sub = paste0('GCM = ', gcm))
plot(envMap, add = TRUE)
plot(rangeMap, border = 'blue', add = TRUE)
map("state", add = TRUE)
map("world", add = TRUE)
points(records$longitude, records$latitude, pch = 16, cex = 0.6, col = 'red')

Fraxinus caroliniana

speciesAb_ <- 'Frax_Caro'
load(paste0('./Models/model_outputs/', speciesAb_, '_GCM', gcm, 
                            '_PC', pc, '.rData'))
plot(rangeMap, border = 'blue', main = substitute(paste('Maxent output, ', 
                                                        italic('Fraxinus caroliniana '),
                                                        'occurrences')),
     sub = paste0('GCM = ', gcm))
plot(envMap, add = TRUE)
plot(rangeMap, border = 'blue', add = TRUE)
map("state", add = TRUE)
map("world", add = TRUE)
points(records$longitude, records$latitude, pch = 16, cex = 0.6, col = 'red')

Fraxinus cuspidata

** Check values in E Texas. plot(raster) extract(raster, click(n=1))

speciesAb_ <- 'Frax_Cusp'
load(paste0('./Models/model_outputs/', speciesAb_, '_GCM', gcm, 
                            '_PC', pc, '.rData'))
plot(rangeMap, border = 'blue', main = substitute(paste('Maxent output, ', 
                                                        italic('Fraxinus cuspidata '),
                                                        'occurrences')),
     sub = paste0('GCM = ', gcm))
plot(envMap, add = TRUE)
plot(rangeMap, border = 'blue', add = TRUE)
map("state", add = TRUE)
map("world", add = TRUE)
points(records$longitude, records$latitude, pch = 16, cex = 0.6, col = 'red')

Fraxinus greggii

speciesAb_ <- 'Frax_Greg'
load(paste0('./Models/model_outputs/', speciesAb_, '_GCM', gcm, 
                            '_PC', pc, '.rData'))
plot(rangeMap, border = 'blue', main = substitute(paste('Maxent output, ', 
                                                        italic('Fraxinus greggii '),
                                                        'occurrences')),
     sub = paste0('GCM = ', gcm))
plot(envMap, add = TRUE)
plot(rangeMap, border = 'blue', add = TRUE)
map("state", add = TRUE)
map("world", add = TRUE)
points(records$longitude, records$latitude, pch = 16, cex = 0.6, col = 'red')

Fraxinus latifolia

speciesAb_ <- 'Frax_Lati'
load(paste0('./Models/model_outputs/', speciesAb_, '_GCM', gcm, 
                            '_PC', pc, '.rData'))
plot(rangeMap, border = 'blue', main = substitute(paste('Maxent output, ', 
                                                        italic('Fraxinus latifolia '),
                                                        'occurrences')),
     sub = paste0('GCM = ', gcm))
plot(envMap, add = TRUE)
plot(rangeMap, border = 'blue', add = TRUE)
map("state", add = TRUE)
map("world", add = TRUE)
points(records$longitude, records$latitude, pch = 16, cex = 0.6, col = 'red')

Fraxinus nigra

speciesAb_ <- 'Frax_Nigr'
load(paste0('./Models/model_outputs/', speciesAb_, '_GCM', gcm, 
                            '_PC', pc, '.rData'))
plot(rangeMap, border = 'blue', main = substitute(paste('Maxent output, ', 
                                                        italic('Fraxinus nigra '),
                                                        'occurrences')),
     sub = paste0('GCM = ', gcm))
plot(envMap, add = TRUE)
plot(rangeMap, border = 'blue', add = TRUE)
map("state", add = TRUE)
map("world", add = TRUE)
points(records$longitude, records$latitude, pch = 16, cex = 0.6, col = 'red')

Fraxinus pennsylvanica

speciesAb_ <- 'Frax_Penn'
load(paste0('./Models/model_outputs/', speciesAb_, '_GCM', gcm, 
                            '_PC', pc, '.rData'))
plot(rangeMap, border = 'blue', main = substitute(paste('Maxent output, ', 
                                                        italic('Fraxinus pennsylvanica '),
                                                        'occurrences')),
     sub = paste0('GCM = ', gcm))
plot(envMap, add = TRUE)
plot(rangeMap, border = 'blue', add = TRUE)
map("state", add = TRUE)
map("world", add = TRUE)
points(records$longitude, records$latitude, pch = 16, cex = 0.6, col = 'red')

Fraxinus profunda

speciesAb_ <- 'Frax_Prof'
load(paste0('./Models/model_outputs/', speciesAb_, '_GCM', gcm, 
                            '_PC', pc, '.rData'))
plot(rangeMap, border = 'blue', main = substitute(paste('Maxent output, ', 
                                                        italic('Fraxinus profunda '),
                                                        'occurrences')),
     sub = paste0('GCM = ', gcm))
plot(envMap, add = TRUE)
plot(rangeMap, border = 'blue', add = TRUE)
map("state", add = TRUE)
map("world", add = TRUE)
points(records$longitude, records$latitude, pch = 16, cex = 0.6, col = 'red')

Fraxinus quadrangulata

speciesAb_ <- 'Frax_Quad'
load(paste0('./Models/model_outputs/', speciesAb_, '_GCM', gcm, 
                            '_PC', pc, '.rData'))
plot(rangeMap, border = 'blue', main = substitute(paste('Maxent output, ', 
                                                        italic('Fraxinus quadrangulata '),
                                                        'occurrences')),
     sub = paste0('GCM = ', gcm))
plot(envMap, add = TRUE)
plot(rangeMap, border = 'blue', add = TRUE)
map("state", add = TRUE)
map("world", add = TRUE)
points(records$longitude, records$latitude, pch = 16, cex = 0.6, col = 'red')

Now, for Lorenz ccsm gcm.

Fraxinus americana

speciesAb_ <- 'Frax_Amer'
load(paste0('./Models/model_outputs/', speciesAb_, '_GCM', gcm, 
                            '_PC', pc, '.rData'))
plot(rangeMap, border = 'blue', main = substitute(paste('Maxent output, ', 
                                                        italic('Fraxinus americana '),
                                                        'occurrences')),
     sub = paste0('GCM = ', gcm))
plot(envMap, add = TRUE)
plot(rangeMap, border = 'blue', add = TRUE)
map("state", add = TRUE)
map("world", add = TRUE)
points(records$longitude, records$latitude, pch = 16, cex = 0.6, col = 'red')

Fraxinus anomala

speciesAb_ <- 'Frax_Anom'
load(paste0('./Models/model_outputs/', speciesAb_, '_GCM', gcm, 
                            '_PC', pc, '.rData'))
plot(rangeMap, border = 'blue', main = substitute(paste('Maxent output, ', 
                                                        italic('Fraxinus anomala '),
                                                        'occurrences')),
     sub = paste0('GCM = ', gcm))
plot(envMap, add = TRUE)
plot(rangeMap, border = 'blue', add = TRUE)
map("state", add = TRUE)
map("world", add = TRUE)
points(records$longitude, records$latitude, pch = 16, cex = 0.6, col = 'red')

Fraxinus caroliniana

speciesAb_ <- 'Frax_Caro'
load(paste0('./Models/model_outputs/', speciesAb_, '_GCM', gcm, 
                            '_PC', pc, '.rData'))
plot(rangeMap, border = 'blue', main = substitute(paste('Maxent output, ', 
                                                        italic('Fraxinus caroliniana '),
                                                        'occurrences')),
     sub = paste0('GCM = ', gcm))
plot(envMap, add = TRUE)
plot(rangeMap, border = 'blue', add = TRUE)
map("state", add = TRUE)
map("world", add = TRUE)
points(records$longitude, records$latitude, pch = 16, cex = 0.6, col = 'red')

Fraxinus cuspidata

** Check values in E Texas. plot(raster) extract(raster, click(n=1))

speciesAb_ <- 'Frax_Cusp'
load(paste0('./Models/model_outputs/', speciesAb_, '_GCM', gcm, 
                            '_PC', pc, '.rData'))
plot(rangeMap, border = 'blue', main = substitute(paste('Maxent output, ', 
                                                        italic('Fraxinus cuspidata '),
                                                        'occurrences')),
     sub = paste0('GCM = ', gcm))
plot(envMap, add = TRUE)
plot(rangeMap, border = 'blue', add = TRUE)
map("state", add = TRUE)
map("world", add = TRUE)
points(records$longitude, records$latitude, pch = 16, cex = 0.6, col = 'red')

Fraxinus greggii

speciesAb_ <- 'Frax_Greg'
load(paste0('./Models/model_outputs/', speciesAb_, '_GCM', gcm, 
                            '_PC', pc, '.rData'))
plot(rangeMap, border = 'blue', main = substitute(paste('Maxent output, ', 
                                                        italic('Fraxinus greggii '),
                                                        'occurrences')),
     sub = paste0('GCM = ', gcm))
plot(envMap, add = TRUE)
plot(rangeMap, border = 'blue', add = TRUE)
map("state", add = TRUE)
map("world", add = TRUE)
points(records$longitude, records$latitude, pch = 16, cex = 0.6, col = 'red')

Fraxinus latifolia

speciesAb_ <- 'Frax_Lati'
load(paste0('./Models/model_outputs/', speciesAb_, '_GCM', gcm, 
                            '_PC', pc, '.rData'))
plot(rangeMap, border = 'blue', main = substitute(paste('Maxent output, ', 
                                                        italic('Fraxinus latifolia '),
                                                        'occurrences')),
     sub = paste0('GCM = ', gcm))
plot(envMap, add = TRUE)
plot(rangeMap, border = 'blue', add = TRUE)
map("state", add = TRUE)
map("world", add = TRUE)
points(records$longitude, records$latitude, pch = 16, cex = 0.6, col = 'red')

Fraxinus nigra

speciesAb_ <- 'Frax_Nigr'
load(paste0('./Models/model_outputs/', speciesAb_, '_GCM', gcm, 
                            '_PC', pc, '.rData'))
plot(rangeMap, border = 'blue', main = substitute(paste('Maxent output, ', 
                                                        italic('Fraxinus nigra '),
                                                        'occurrences')),
     sub = paste0('GCM = ', gcm))
plot(envMap, add = TRUE)
plot(rangeMap, border = 'blue', add = TRUE)
map("state", add = TRUE)
map("world", add = TRUE)
points(records$longitude, records$latitude, pch = 16, cex = 0.6, col = 'red')

Fraxinus pennsylvanica

speciesAb_ <- 'Frax_Penn'
load(paste0('./Models/model_outputs/', speciesAb_, '_GCM', gcm, 
                            '_PC', pc, '.rData'))
plot(rangeMap, border = 'blue', main = substitute(paste('Maxent output, ', 
                                                        italic('Fraxinus pennsylvanica '),
                                                        'occurrences')),
     sub = paste0('GCM = ', gcm))
plot(envMap, add = TRUE)
plot(rangeMap, border = 'blue', add = TRUE)
map("state", add = TRUE)
map("world", add = TRUE)
points(records$longitude, records$latitude, pch = 16, cex = 0.6, col = 'red')

Fraxinus profunda

speciesAb_ <- 'Frax_Prof'
load(paste0('./Models/model_outputs/', speciesAb_, '_GCM', gcm, 
                            '_PC', pc, '.rData'))
plot(rangeMap, border = 'blue', main = substitute(paste('Maxent output, ', 
                                                        italic('Fraxinus profunda '),
                                                        'occurrences')),
     sub = paste0('GCM = ', gcm))
plot(envMap, add = TRUE)
plot(rangeMap, border = 'blue', add = TRUE)
map("state", add = TRUE)
map("world", add = TRUE)
points(records$longitude, records$latitude, pch = 16, cex = 0.6, col = 'red')

Fraxinus quadrangulata

speciesAb_ <- 'Frax_Quad'
load(paste0('./Models/model_outputs/', speciesAb_, '_GCM', gcm, 
                            '_PC', pc, '.rData'))
plot(rangeMap, border = 'blue', main = substitute(paste('Maxent output, ', 
                                                        italic('Fraxinus quadrangulata '),
                                                        'occurrences')),
     sub = paste0('GCM = ', gcm))
plot(envMap, add = TRUE)
plot(rangeMap, border = 'blue', add = TRUE)
map("state", add = TRUE)
map("world", add = TRUE)
points(records$longitude, records$latitude, pch = 16, cex = 0.6, col = 'red')